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1.  INTRODUCTION 


Requirements  for  advanced  gun  systems,  such  as  the  U.S.  Anny  Advanced  Field  Artillery 
System  (AFAS),  include  increased  firing  rates  in  both  burst  and  sustained  firing  modes.  A 
serious  problem  attendant  to  such  capabilities  is  that  of  thermal  management,  which  is 
required  to  minimize  thermal  damage  and  wear  to  the  tube  and  associated  gun  parts,  to 
reduce  the  effects  of  thermai  distortions  on  gun  tube  accuracy,  and  to  eliminate  safety 
problems  such  as  charge  or  projectile  cookoff.  A  comprehensive  approach  to  the  problem 
Involves  a  careful  analysis  of  the  influence  of  firing  rate  on  both  the  weapon  heating  problem 
and  mission  effectiveness  of  each  of  tiie  propulsion  options  under  consideration.  In  the  case 
of  AFAS.  this  choice  now  includes  solid  and  liquid  propellants. 

A  two*dimensional  (2-D)  gun  tube  heat  conduction  code  was  developed  especially  to 
address  the  current  thermal  msuiagement  concerns  (Chandra  and  Rsher  1989a.  1989b; 
Chandra  1990;  Conroy  1990;  Keller  et  al.  1991).  Tills  code  uses  the  gas  tempeiature, 
pressure,  and  velocity  output  of  an  interior  ballistics  (IB)  code  as  its  input.  It  is  assumed  that 
the  heat  transfer  process  can  be  uncoupled  from  the  combustion  and  flow  occurring  within  the 
gun  tube.  The  interior  ballistics  code  provides  a  file  containing  gas  temperature,  gas 
pressure,  and  gas  veioci^  as  a  function  of  time  and  distance  along  the  gun  tube.  The  IB 
calculation  is  continued  until  the  gas  temperature  and  pressure  dedine  to  near  ambient  (i.e., 
the  blowdovm  phase  must  be  included). 

For  a  solid  propellant  charge,  ttie  XNCVAKTC  code  (Gough  1980, 1990)  is  used  to 
generate  the  required  profiles.  This  report  discusses  the  code  developed  to  run  the  liquid 
propellant  simulations. 

2.  LIQUID  GUN  COOES 

Most  of  the  regenerative  liquid  propellant  gun  modeling  at  the  U.S.  Army  Ballistic 
Research  Laboratory  (BRL),  Aberdeen  Proving  Ground,  MD,  has  used  a  lumped  parameter 
code  developed  in-house  (Coffee  1985, 1988;  Mom'son  and  Coffee  1990).  This  code  has 
been  extensively  validated  agednst  experimental  data  (Coffee.  Wren,  and  Morrison  1 989, 
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1990;  Wren,  Coffee,  and  Morrison  1990;  Coffee  et  al.  1991).  Some  special  purpose  codes 
have  also  been  developed— a  one-dimensional  (l-D)  modsi'  for  the  liquid  reservoir 
(Coffee  1988)  and  a  2-D  model  for  the  combustion  chamber/gun  tube  (Coffee,  to  be 
published).  Also,  Gough  has  written  several  liquid  propellant  models— a  model  that  is  1-0  in 
the  gun  tube  and  lumped  parameter  elsewhere  (Gough  1987),  and  a  complete  1-0  model 
(Wren  and  Gough  1991).  However,  the  solution  terminates  at  projectile  exit. 

The  new  factor  required  is  the  modeling  of  the  blowdown  phase.  After  the  projectile  exits 
the  gun  tube,  a  rarefaction  wave  will  propagate  from  the  muzzle  toward  the  breech.  It  is 
possible  to  model  a  rarefaction  wave  using  a  lumped  parameter  code.  In  fact,  the  lumped 
parameter  code  assumes  that  a  raretecti  jn  wave  propagates  from  the  chamber  toward  the 
muzzle  after  bum  out  of  the  propellant  (Morrison  and  Coffee  1990).  However,  rather  than 
modeling  two  rarefaction  waves  with  a  lumped  parameter  code,  a  1-D  model  In  the  gun  tube 
wsfi  enveloped.  For  the  rest  of  the  system,  the  lumped  parameter  representation  Is 
sufficlentiy  accurate. 

3.  NUMERICAL  PROCEDURE 

In  the  Concept  VIC  regenerative  liquid  gun,  the  propellant  is  initially  behind  two  pistons. 
The  inner  or  control  piston  moves  first,  opening  up  an  annular  vent  The  outer,  or  injection, 
piston  then  slightly  trails  the  inner  piston.  The  motion  of  the  control  piston  is  modulated  by  a 
damper  region.  Propellant  is  injected  between  the  pistons  Into  the  combustion  chamber.  The 
combustion  takes  place  in  the  combustion  chamber,  and  the  gas  then  flows  into  the  gun  tube. 
Unlike  solid  propellant  guns,  there  is  a  large  area  change  from  the  chamber  to  the  tube. 

The  governing  equations  from  the  lumped  parameter  model  are  used  for  the  damper, 
liquid  reservoir,  combustion  chamber,  and  projectile.  Previous  work  indicates  that  very  little 
liquid  gets  into  the  gun  tube,  so  for  simplicity,  the  fluid  in  the  tube  is  assumed  to  be  pure  gas. 

A  staggered  grid  Is  used  in  the  gun  tube  (Figure  1).  The  velocities  are  defined  at  the  solid 
lines.  This  vector  grid  Is  evenly  spaced  between  the  gun  tube  entrance  and  the  projectile. 

The  scalar  quantities  (pressure,  temperature,  and  density)  are  defined  midway  between  the 
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velocities  (dotted  lines).  The  grid  divides  the  tube  into  scalar  control  volumes  (between  the 
solid  lines).  The  scalar  quantities  are  assumed  to  be  uniform  in  the  scalar  control  volumes, 
with  a  discontinuity  at  the  solid  lines.  Similarly,  tiiere  are  vector  control  volumes  (between  the 
dotted  lines).  The  velocity  is  assumed  to  be  uniform  within  vector  control  volumes.  Additional 
control  volumes  are  established  outside  the  physical  gun  tube  to  help  set  the  boundary 
conditions  (see  below).  The  projectile  motion  is  computed  in  the  usual  fashion,  based  on  the 
pressure  at  the  base  of  the  projectile.  As  the  projectile  moves,  the  grid  stretches  (conformal 
gild).  Additional  ooints  are  added  to  the  grid  as  needed. 


The  governing  equations  are  cast  in  integral  form.  The  procedure  is  inherently 
conservative.  That  is,  if  mass,  momentum,  or  energy  leaves  one  control  volume  through  a 
boundary,  it  automatically  enters  the  neighboring  control  volume.  For  simplicity,  an  explicit 
numerica'  integration  scheme  is  used.  Upwind  differencing  is  used  for  stability. 

3.1  Governing  Equations.  The  1-D  continuity  equation  is  (Bird,  Stewart,  and 
Ughtfoot  1960) 


9p/3f  +  3(p  v)/3x  ■  0  , 


(1) 


where  p  Is  the  gas  density  (g/cm^  and  v  is  the  gas  velocity  (cm/s).  Now  consider  some 
volume  In  three-dimensional  (3-D)  space.  By  the  divergence  theorem.  Equation  1  can  be 
written  as 


3/3fJ^cf\^-- J^(3(pv)/3x)  dV 

•  - n  dS •  (2) 

The  first  two  integrals  are  volume  integrals.  The  last  integral  is  a  surface  integral.  The  vector 
I  Is  the  unit  vector  in  the  x  (down  bore)  direction,  and  n  Is  the  unit  outward  normal  vector  to 
the  surface.  So  any  change  in  the  mass  in  the  control  volume  is  due  to  mass  flux  across  a 
boundary. 

The  energy  equation  can  be  written  as  (Bird,  Stewart,  and  Ughtfoot  1960) 


p  c,  3  r/3f  -  - p  c,  V  3 T/3x  -  T  {dp/d T)^  ( 3 v/3x) ,  (3) 


where  Tis  the  temperature  (K),  is  the  specific  heat  at  constant  volume  (j/g  -  K),  and  p  is 
the  pressure  (MPa).  The  viscous  heating  term,  which  is  almost  always  small,  is  eliminated. 
The  thermal  conductivity  term  is  also  not  included.  In  the  gun  tube,  the  temperature  is  fairly 
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constant,  and  the  thermal  conductivity  term  is  much  smaller  than  the  convection  terms.  Using 
the  Noble-Abel  equation,  the  derivative  of  p  with  respect  to  T  with  tie  density  constant  is 
simply  the  pressure  over  the  temperature.  Making  this  substitution  and  dividing  by  the  specific 
heat  produces 


p  9779f  ■ -p  /  97V9x  -  p/Cj,  (9v/9x).  (4) 

Using  the  continuity  equation,  this  can  be  written  as 


9(p  T)/9t- -9{pvr)/9x-p/c,(9vr/9x).  (5) 


The  divergence  theorem  Is  applied  The  pressure  is  assumed  to  be  constant  in  the  control 
volume  and  can  be  taken  outside  the  integral.  Then 


9/9fJ^(pT)dV-  - J^iQvTI)n  dS-  (p/c,)  dS  (6) 


The  basic  variable  is  the  mass  times  the  temperature.  This  can.  as  before,  be  changed  by 
convection  (first  right-hand  term).  Hot  or  cold  gas  can  enter  through  a  boundary.  The 
temperature  can  also  be  changed  through  work  (second  right-hand  term).  Note  that  in  the  last 
term,  the  pressure  is  not  evaluated  at  the  boundary,  but  is  the  pressure  in  the  control  volume. 

The  momentum  equation  in  the  axial  direction  is  (Bird,  Stewart,  and  Ughtfoot  1960) 


p  dv/dt  •  -pv  9v/9x  -  9p/9x. 


(7) 


The  viscosity  terms  are  not  included.  Viscosity  will  only  be  important  for  shear  stress 
(perpendicular  to  the  direction  of  flow).  For  the  present  1-0  model,  this  cannot  be  resolved. 
Using  the  continuity  equation,  this  can  be  rewritten  as 
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d{pv)/dt  ■  -d(pv^)/dx  -  dpIdX. 


(8) 


Applying  the  divergence  theorem, 


dfdtj^{pv)clV^  I  +  pf)  n  dS. 


(9) 


32  Grid.  Figure  1  shows  the  grid  used  in  the  gun  tube.  The  solid  lines  represent  the 
vector  grid  and  are  labeled  from  1  to  nt  The  grid  points  are  evenly  spaced.  Midway  between 
the  vector  grid  u  the  scalar  grid,  labeled  from  1  to  nt  - 1. 

The  vector  grid  splits  the  chamber  up  Into  scale:  control  volumes.  The  scalar  quantities 
are  assumed  to  be  constant  vrithin  each  scalar  volume  (between  solid  lines).  Ibere  is  a 
discontinuity  bt  the  scalar  quantities  at  each  solid  line.  Each  scalar  volume  has  a  volume  V 
and  a  gas  mass  M  The  gas  in  the  scalar  volume  hfti  density  p  «  MV.  pressure  p,  and 
temperature  T, 

The  dotted  lines  spUt  up  the  chamber  btto  vector  control  volumes.  The  velocities  are 
assumed  to  be  the  same  through  each  vector  control  volume,  and  there  is  a  jump  in  velocity 
at  the  dotted  lines.  Each  vector  control  volume  Is  made  up  of  half  of  two  scalar  contrcM 
volumes.  The  mass  in  the  vector  control  volume  /  is  hsM  the  mass  in  the  scalar  control 
volume  /  -  f  plus  half  the  mass  in  the  scalar  control  volume  /,  that  is,  M^l)  «  0.5  (M(l  ~1) 

•f  M(0).  Ute  vector  contrr^  volumes  at  the  boundary  require  special  consideration  (see 
Section  3.4).  If  the  grid  spacing  exceeds  a  user-spedtied  value,  dmax,  an  additional  grid 
point  is  added.  A  maximum  number  of  grid  points,  nmax,  can  also  be  specified. 

The  grid  stretches  between  the  tube  entrance  and  the  projectile.  The  grid  velocity  U(nt)  at 
the  projectile  equals  the  projectile  veiodty.  The  grid  veiodty  U(1)  at  the  tube  entrance  Is  zero. 
The  other  grid  velocities  are  found  by  linear  interpolation.  So  the  grid  remains  equally  spaced 
as  the  Integration  proceeds.  At  muzzle  exit  the  grid  is  attached  to  the  end  of  the  gun  tube 
and  all  the  grid  velocities  are  set  equal  to  zero. 
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With  this  formulation,  the  governing  equations  can  be  simplified.  Consider  the  continuity 
Equation  2  applied  over  a  scalar  control  volume.  The  density  is  constant  within  the  control 
volume,  so  It  can  be  taken  outside  of  the  left-hand  integral.  The  left-hand  side  then  just 
becomes  density  times  volume  or  mass.  The  right-hand  side  involves  the  velocity  across  the 
boundaries.  Since  the  boundaries  are  also  moving,  this  is  the  difference  between  the  grid  and 
gas  velocities.  The  density  must  also  be  known  on  the  boundaries.  Upwind  differencing  is 
used,  which  means  the  density  on  the  boundary  is  assumed  to  be  the  density  from  the 
direction  of  flow.  The  direction  of  flow  on  a  boundary  is  determined  by  the  sign  of  the  gas 
minus  grid  velocity,  defined  as 


u{l)  -  v{l)  -  U{l). 


(10) 


Assuming  a  positive  gas  minus  grid  velocity  and  using  upwind  differencing, 


aW(/)/8f-{p(/-1)  w(/) -p(/)  y(/+1)}A  0^) 


where  A  Is  the  gun  tube  area  (cm^.  If  instead  the  quantity  u(l)  is  negative,  the  corresponding 
boundary  densities  are  p^i;  and  p(l+  1). 

Similariy,  the  energy  equation  (6)  becomes 


a[Af7](/)/af-{p{/- 1)  r(/- 1)  u(/)  -p(/)  r(/)u(/  + 1))  a 

MP(/)/<v1M/)- WZ+I)}  A  . 


(12) 


The  first  term  on  the  right-hand  side  is  a  convection  term,  and  the  appropriate  velocity  is  the 
velocity  across  the  boundary.  The  second  term  Is  a  work  term.  The  change  in  temperature  is 
due  to  the  compression  or  expansion  of  the  gas  and  the  appropriate  velocity  is  the  gas 
velocity. 
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d[M,v]{i)m  «  {p(/  -  1)  v(/  -  1)  u{i  -  1)  -  p(/)  v{i)  u{i))  A 
*9o  <P(/  -  1)  -PiO^  A  , 


(13) 


where  Is  a  conversion  constant  (10'  gm/s-cm-MPa).  For  this  case,  the  boundaries  of  a 
vector  control  volume  are  in  tiie  middle  of  ttie  scalar  Cjntrol  volumes,  so  the  density  and 
pressure  are  ttie  same*  Irregardless  of  the  direction  of  flow.  However,  tiie  velocity  must  be 
found  on  the  veiodty  control  volume  boundary,  suid  s^jain  upwind  differencing  is  used.  The 
form  given  In  Equation  13  Is  for  flow  In  the  positive  direction,  whidi  will  be  the  case  until  late 
in  the  blowdown  phase. 

3.3  Heat  Loss.  The  loss  of  energy  to  the  gun  tube  walls  has  a  noticeable  effect  on  tiie 
interior  ballistics.  To  compute  heat  transfer  to  the  wail,  a  2-D  code  is  required,  instead,  a 
correlation  from  the  tumped  parameter  code  is  implemented. 

The  gun  tube  temperature  Is  assumed  to  remain  const^  (infinite  sink).  Then  the  heat 
loss  can  be  represented  as 


Q^m4h„{T-TJ/d  (14) 

where  T  is  the  gas  temperature  In  the  tube,  d  Is  the  tube  dameter  (cm),  Is  the  heat  transfer 

coefficient  (j/cm^  <•  K  s),  and  Q,  is  the  heat  loss  (j/cm’  -  s)  (Bird,  Stewart  and  Ughtfoot 
I960).  The  variation  of  teir^rature  in  the  gun  tube  is  ignored.  The  dfficutty  is  in  computing 
the  heat  loss  coefficient  in  tiie  code,  ffie  correlation  of  Sleder  and  Tate  for  pipe  flow  is  used. 
This  assumes  hilly  developed  flow  in  a  smooth  pipe  with  a  nearly  constant  wall  temperature. 
For  high^  turbulent  flow. 


-  (Wd)  0.026  ,  (15) 


where  k\s  the  tiiermai  conductivity  of  the  gas  (j/cm  -  K  -  s).  Re  is  the  Reynold’s  number, 
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Pr  is  the  Prandtl  number,  and  v  Is  the  viscosity  (poise  =  g/cm  -  s).  The  correlation  is  for 
steady-state  flow,  while  the  gun  conditions  are  highly  transient  More  importantly,  the 
correlation  has  been  found  to  be  accurate  for  Reynold's  numbers  up  to  about  10’,  while 
Reynold's  numbers  in  the  gun  tube  can  be  two  orders  of  magnitude  higher. 

if  a  HAN-based  propellant  is  completely  burned,  it  will  form  a  mixture  of  COj,  HjO,  and  N^. 
Procedures  have  been  developed  to  find  ttie  viscosity  and  thermal  conductivity  of  a  mixture  of 
gases  (Coffee  and  Helmeri  1981).  The  procedures  were  developed  for  tiie  low-density  limit 
Fortunately,  for  high  temperatures,  the  low-density  limit  is  veiy  accurate  even  for  very  high 
pressures  (Bird,  Stewart  and  Ughtfoot  1960).  Rts  were  made  for  the  viscosity  and  thermal 
conductivity  cf  the  gases  resulting  from  HAN1845  between  300  K  and  3,000  K.  Values  for 
other  HAN-based  propellants  are  almost  tiie  same.  A  separate  low  temperature  fit  was  made 
to  detamfiine  the  viscosity  at  the  gun  tube  wall  The  fits  are 


v^3275  10-« 
V,  •  5.3083  10’*^ 

it  -  1.7974  10-* 


(16) 


The  Reynold's  number  is  given  by 


Re  -  p  vdV 


(17) 


and  the  Prandtl  number  by 


Pr  -  Cp  v//r  , 


(18) 


where  Cp  Is  the  spedfic  heat  at  constant  pressure. 
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In  the  lumped  parameter  code,  the  correlation  Is  applied  to  the  tube  as  a  whole,  using 
space  mean  values  for  pvand  T.  In  the  1-D  code,  the  correlation  Is  applied  to  each  scalar 
control  volume. 

Note  that  In  ttie  actual  thermal  management  code,  a  different  correlation  is  used  for  heat 
transfer,  and  the  gun  tube  can  change  temperature.  The  effects  of  these  inconsistencies  has 
not  yet  been  studied. 

3.4  Boundary  Conditions.  As  usual,  the  ftow  into  the  gun  tube  is  assumed  to  be 
steady-state  isentropic  flow  (Coffee  1988).  This  gives  the  mass  flow  Into  the  tube.  A  scalar 
control  volume  Is  set  up  just  to  the  left  of  the  gun  tube  (I «  0).  The  density  and  pressure  In 
this  control  volume  is  set  equal  to  the  density  and  pressure  Just  inside  the  tube  (at  scalar 
control  volume  I  •  1).  The  temperature  is  set  equal  to  the  chamber  gas  temperature,  so  that 
mass  times  temperature  is  conserved  by  flow  across  the  boundary.  Then  the  gas  velocity  at 
the  entrance  is  given  by 


v(1).mf/[p(0)A]  , 


(19) 


where  m/is  the  mass  flux  (g/s). 

Before  muzzle  exit,  the  velocity  at  the  right-hand  side  of  the  tdse  Is  Just  the  projectile 
velocity.  For  the  scalars,  a  reflection  condition  is  imposed  (i.e.,  the  densify,  temperature,  and 
pressure  Just  outside  the  region  of  integration  0  «  nt]  Is  equal  to  the  scalars  Just  at  the  end  of 
the  physical  tube  p  a  nt  ~  1  ]). 

i 

I 

At  muzzle  exit,  the  grid  is  attached  to  the  end  of  the  gun  tube.  The  flow  out  of  the  tube  is 
assumed  to  be  steac^-stats  isentropic  flow.  If  the  outflow  velocity  is  tevger  than  the  local 
sound  speed.  It  must  be  lowered  to  the  sonic  velocity  (choked  flow).  The  flow  will  be  choked 
for  most  of  the  blowdown  phase. 

For  isentropic  flow,  the  relevant  equations  are  the  following.  The  first  law  of 
thermodynamics 
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h  +  v^/2g„^  cons^t , 


(20) 


where  h  Is  the  gas  enthalpy  (j/g)  and  ir  Is  the  gas  velocity  (cm/s).  The  Isentropic  process 
equation 


p  (1  /p  -  by  •  constant , 

where  p  Is  the  pressure  (MPa),  p  is  the  density  (g/om^,  d  is  the  oovoiume  (cmVg),  and  y  Is 
the  ratio  of  specific  heats.  The  equation  of  state  is 

pH,  T.p(1 -6p),  (22) 


where  7is  the  temperature  (K),  and  is  the  specific  gas  constant  (j/g  -  K).  Aiso  required  is 
the  equation  for  the  speed  of  sound  in  a  gas, 

e*-fcYH/tP(l -ip)]  .  (2®) 


and  the  enthaipy  equation. 


h»c^T*  bp,  (24) 

where  c^ls  the  specific  heat  at  constant  pressure  (j/g  -  K). 

Let  p, ,  p, ,  and  r,  be  the  scalar  quantities  in  the  scaiar  controi  volume  at  the  end  of  the 
gun  tube,  and  v,  the  gas  veiodty  at  the  vector  grid  point  just  before  the  muzzle  exit  The 
problem  Is  then  finding  the  quantities  p, ,  p, ,  7, ,  and  v,  at  the  throat  The  velocity  v,  is 
assumed  to  be  the  velocity  exactly  at  the  end  of  the  tube.  The  scalsLr  quantities  are  tfie 
values  In  the  control  volume  just  outsi  la  the  tube. 
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Suppose  the  throat  pressure  is  known.  The  throat  density  is  found  by  the  process 
equation.  The  temperature  is  found  from  the  equation  of  state.  Then  from  Equations  20 
and  24, 


(25) 


Previously,  the  interest  has  been  in  subsonic  flow.  In  this  case,  the  throat  pressure  can  just 
be  taken  as  the  pressure  just  outside  the  throat  (in  this  case,  as  1  atm).  However,  when  the 
flow  reaches  the  local  speed  of  sound  (in  the  throat),  the  flow  becomes  choked,  and  no  further 
increase  in  the  pressure  drop  will  have  any  effect  The  gas  cannot  move  away  fast  enough, 
and  a  higher  pressure  region  will  exist  just  outside  the  gun  tube.  The  pressure  at  which  the 
flow  becomes  choked  Is  called  the  critica!  pressure. 

For  an  Ideal  gas.  the  critical  pressure  can  be  found  analytically.  Let  the  mach  number 
M  >  v/a  Define  the  stagnation  pressure  p.  as  the  pressure  obtained  when  a  gas  is 
decelerated  Isentropically  to  zero  velocity.  Then  (Fox  and  McDonald  1985) 

p.  -  p  (1  +  (y  -  1 )  (26) 


Tlten  the  critical  pressure  p^  is  given  by 


(27) 


The  procedure  is  to  And  the  critical  pressure  based  on  the  conditions  at  the  right  end  of  the 
tube.  If  the  critical  pressure  is  less  than  the  exit  pressure,  the  flow  is  subsonic,  and  the  exit 
pressure  is  used  as  the  throat  pressure,  if  the  criticai  pressure  is  greater  than  the  exit 
pressure,  the  flow  is  choked,  and  the  criticai  pressure  is  used  as  the  throat  pressure. 

For  the  Noble-Abel  equation  of  state,  the  relevant  equation  becomes  transcendental,  and 
the  critical  pressure  cannot  be  found  analytically.  Instead,  an  iterative  procedure  is  used.  The 
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ideal  gas  law  critical  pressure  is  computed  as  a  starting  point.  Using  this  as  the  throat 
pressure,  p, ,  p, ,  T, ,  v,,  and  c,  are  computed.  If  the  gas  velocity  v^  is  less  than  the  sound 
speed,  pf  is  decreased  by  an  amount  proportional  to  the  difference  in  the  velocities.  If  v,  is 
greater  than  c, ,  p,  is  increased.  The  ideal  gas  law  critical  pressure  is  a  good  approximation, 
so  convergence  is  rapid. 

There  is  one  last  problem.  The  mass  flow  out  of  the  gun  tube  should  be  p,  v,  A.  where  A 
is  the  area  of  the  tube.  However,  the  code  uses  upwind  differencing,  so  the  normal  procedure 
would  give  a  mass  flow  of  p,  v,  A.  So  special  logic  is  included  at  the  end  of  the  tube  to  obtain 
the  proper  mass  flux  out. 

Late  in  the  blowdown  phase,  the  pressure  inside  the  tube  can  become  lower  than 
atmospheric  pressure.  In  that  case,  there  is  flow  into  the  tube.  The  gas  outside  the  tube  is 
assumed  to  be  stagnant,  since  there  is  no  simple  way  to  compute  its  velocity. 

On  the  other  hand,  flow  from  the  chamber  into  the  tube  is  assumed  to  always  be  positive, 
since  back  flow  has  not  previously  been  important,  if  the  conditions  indicate  back  flow,  the 
mass  flux  is  just  set  to  zero.  This  has  not  been  considered  important  enough  to  update. 

3.5  Time  Step.  To  take  a  time  step,  the  time  derivatives  of  the  lumped  parameter 
quantities  are  firjit  evaluated  and  stored  in  a  vector  YOOT.  The  time  step  ctf  is  a  user- 
specified  parameter.  Copies  of  the  lumped  parameter  values  are  stored  in  a  vector  Y.  This  is 
updated  using  the  Euler  formula. 


Y  -  Y  +  YDOT*dt  .  (28) 

However,  the  actual  lumped  parameter  values  are  not  updated,  so  the  code  Is  fully  explidt. 

To  take  a  tube  time  step,  the  boundary  conditions  are  applied  to  find  the  velocity  at  the 
tube  entrance  and  at  the  projectile.  Next,  the  heat  loss  term  is  applied  to  the  energy  equation, 
i.e., 
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(29) 


MT=MT  -  Q^Vdt/c,  . 

Other  than  heat  loss  to  the  walls,  thermal  conductivity  is  neglected. 

Then  the  Lagrange  part  of  the  time  step  is  taken.  That  is,  the  grid  is  momentarily 
assumed  to  move  at  the  same  velocity  as  the  gas.  The  continuity  and  energy  equations  are 
unchanged,  since  there  is  no  mass  flow  through  the  boundaries.  For  the  momentum 
equation,  the  work  term  is  computed 

lM,v]^l)~g,lp^l-^)-p(,l))  A.  (30) 

The  mass  M,  In  the  vector  volume  remains  unchanged.  The  velocities  are  updated  by 

V{l)  (31) 

The  boundary  conditions  are  then  updated. 

The  scalar  control  volumes  have  a  new  Lagrange  volume, 

V, (/).!/(/) +{v(/  + 1)  -  v(/)}  4  df.  (32) 

The  new  density  is  given  by 

p{l)  M  M{i)/V,{l).  (33) 
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The  temperature  equation  is  updated  using  the  work  term 


[MTni).lMTm-ip{l)/c,]lV,-V]  ,  (34) 


and  the  new  temperature  is  now 

7(/)  -  [M  T]{l)/M(l)  .  (35) 

Finally,  the  grid  is  moved  to  Its  proper  position.  The  values  needed  on  the  boundaries  can 
now  be  found  using  upstream  differencing.  Assuming  that  u(l)  is  positive  everywhere,  the 
advection  part  of  the  equations  are 


M{l)  -  M(t)  +  {p{/-1 )  u{l)  -  p(/)  t/(/+1)i  A  (36) 

[WT]  (f) -[M  71(0 +(p  (/-I)  r(M)  a(/) -p{/)  T(/)  u(/  + 1)}  A  (37) 

IH  v']  (/)  -  IH  v'l  (/)  +  <P(/  -  1)  ‘'(Z  -  1)  w(/  -  1)  -  p(0  v(0  u{l))  A.  (38) 

The  boundary  conditions  are  called  again  and  the  scalars  are  updated.  The  pressure  is 
obtained  from  the  equation  of  state. 

Next,  the  grid  is  updated.  The  new  position  of  the  grid  points  is  found  by  applying  the  grid 
velocity  over  the  time  step.  The  size  of  the  new  control  volumes  is  computed.  The 
acceleration  on  the  projectile  is  found  from  the  pressure  at  the  right-hand  scalar  control 
volume.  From  the  acceleration,  the  new  projectile  velocity  is  found.  This  then  allows  the  new 
grid  velocities  to  be  computed.  If  the  control  volumes  are  longer  than  a  user-specified  length, 
dxmax,  a  new  grid  point  is  added.  The  grid  points  are  rearranged  so  the  grid  spacing  is 
unifonn.  The  velocities  on  the  new  grid  are  found  by  interpolation.  The  new  scalar  control 
volumes  are  each  made  up  of  parts  of  two  old  scalar  control  volumes.  The  appropriate 
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fraction  of  material  from  each  of  the  old  control  volumes  is  put  into  the  new  control  volumes, 
and  the  boundary  conditions  are  updated.  Finally,  the  new  values  of  the  lumped  parameter 
variables  are  updated. 


There  is  a  minor  numerical  problem  if  the  projectile  initially  has  zero  offset.  In  the  lumped 
parameter  code,  if  the  projectile  Is  only  a  short  distance  downtube,  the  chamber/tube  Is 
treated  as  one  region.  This  Is  no  longer  desirable,  since  the  chamber  can  have  liquid  and  the 
tube  cannot  So  instead,  a  minimum  offeet  of  0.01  cm  is  assumed.  However,  suppose  a 
moderate-sized  time  step  is  taken.  The  mass  flux  Into  the  tube  is  constant  over  the  time  step. 
Since  the  tube  volume  is  tiny,  the  density  and  pressure  will  overshoot,  leading  to  numerical 
problems.  So  logic  was  added  to  allow  a  smaller  time  step  if  the  grid  spacing  is  smaller  than 
dxmax. 

4.  VALIDATION 

4.1  Comparison  With  Lumped  Parameter  Model.  An  initial  comparison  was  made  with  the 
lumped  parameter  gun  code  (Table  1).  This  can  only  be  done  up  to  muzzle  exit 


Table  1.  Comparison  Between  Lumped  Parameter  and  1-D  Models  for  Round  77 


dt 

(ms) 

dxmax 

(cm) 

nmax 

Vm 

(mfo) 

Maximum 

Chamber 

Pressure 

(MPa) 

Run 

time 

(mln/s) 

Lumped  parameter 

896 

238 

12 

.01 

5 

SO 

917 

256 

16 

.003 

5 

50 

900 

238 

41 

.001 

5 

50 

900 

238 

1:53 

.003 

2 

200 

906 

238 

157 

.001 

2 

200 

902 

239 

527 

.0003 

2 

200 

901 

238 

17:32 
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As  a  test  case,  a  recent  14-liter  (full  charge)  firing  from  the  General  Electric  second 
generation  155-mm  gun  was  used  (Round  77)  (Coffee  and  Wren,  to  be  published).  This 
results  in  a  relatively  high  muzzle  velocity  (893  m/s),  where  the  differences  in  the  gun  tube 
model  are  expected  to  be  more  important  The  projectile  travel  is  just  under  7  m. 

Several  runs  were  made  to  check  on  the  convergence  properties.  The  numerical  controls 
are  the  time  step  dt,  the  grid  spacing  dxmax,  and  the  maximum  number  of  grid  points  nmax. 
The  calculated  muzzle  velocities  and  maximum  chamber  pressures  are  given.  The 
calculations  were  performed  on  a  Silicon  Graphics  4D/310  QTX  workstation.  The  results  show 
convergence  even  with  a  fairly  coarse  grid. 

Figure  2  compares  the  chamber  pressure  from  the  lumped  parameter  model  and  from  the 
two  grids  used  with  the  1>D  model  (smallest  time  steps  from  Table  1).  The  profiles  are  almost 
overlays.  It  Is  only  at  late  times,  when  the  projectile  is  far  downbore,  that  the  1-0  model 
shows  differences  from  the  lumped  parameter  model.  Rgure  3  shows  the  corresponding 
curves  for  the  projectile  velocities. 

Figures  4-7  compare  the  gun  tube  profiles  near  the  end  of  the  firing  cyde  at  18  ms. 
Overail,  the  two  1-0  models  are  in  very  dose  SKjreement,  indicating  that  grid  independence 
has  been  reached.  The  1>D  code  shows  higher  velocities,  but  the  velodty  at  the  projectile  is 
aimost  the  same  as  from  the  iumped  parameter  modei.  The  pressure  curve  from  the  1-0 
modei  shows  lass  curvature  than  the  lumped  parameter  model.  However,  both  models  predict 
a  pressure  diff^rance  of  approximately  5  MPa  from  the  tube  enhance  to  the  projectile.  The 
density  and  temperature  profiles  in  the  gun  tube  are  somewhat  more  complicated.  There  Is  a 
low-density  region  behind  the  projectiier  since  the  motion  of  the  projectile  expands  the  gas. 
The  pressure  equilibrates  faster  than  the  density,  so  the  temperature  goes  up.  On  the  other 
hand.  In  a  slower  projectile  case,  there  is  instead  a  high-density  region  just  behind  the 
projectile.  The  flow  from  the  chamber  tends  to  increase  the  density,  and  the  motion  of  the 
projectile  decreases  the  density.  The  competition  between  these  two  effects  can  lead  to 
somewhat  complicated  behavior  near  the  projectile. 
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Figure  6.  Rourid  77.  Gun  Tube  Gas  Density.  Time  « 18  ms.  Lumped  Parameter  Model 
(Une).  Coarse  Grid  (Pott.  Rnar  Grid /DasM. 


Rgure  7.  Round  77.  Gun  Tube  Gas  Temperature.  Time  =  18  ms.  Lumped  Parameter  Model 
(Line).  Coarse  Grid  (Dot).  Finer  Grid  (Dash). 
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Comparisons  were  also  made  for  a  30-mm  system,  a  105-mm  system,  and  smaller  charge 
155>mm  systems.  The  agreement  between  the  lumped  parameter  code  and  the  1-D  code 
was  better  than  that  reported  above. 

4,2  Comparison  With  Theory.  Experimental  data  for  the  blowdown  phase  is  not  available 
for  liquid  propellant  guns.  So  a  comparison  Is  made  with  a  theoretical  model  developed  by 
Comer  (1950). 

The  Comer  model  assumes  that  the  propellant  is  totally  burned  before  the  projectile  starts 
moving.  The  gas  then  expands  adiabatically  until  muzzle  exit  It  Is  assumed  that  tiie 
pressure  follows  a  standard  Lagrange  gradient  A  rarefaction  wave  develops  and  moves  back 
toward  the  breech.  Comer  developed  expressions  for  the  profiles  from  the  breech  to  the 
rare^ction  wave,  as  well  as  a  formula  for  ttte  location  of  the  rarefaction  wave.  There  is  not  a 
general  formula  for  the  profiles  from  the  rarefaction  wave  to  ttie  muzzle.  Comer  uses 
expsmsions  in  ttie  quantity  charge  over  projectile  mass,  and  ignores  second-order  terms,  so 
his  expressions  are  ortiy  approximate.  The  model  requires  as  Input  the  projectile  muzzle 
velocity  and  the  average  temperature  of  the  gas  in  the  tube  at  muzzle  exit 

in  order  to  mimic  the  assumptions  of  the  Comer  model  In  the  1-0  code,  the  mass  flow 
from  the  chamber  to  the  tube  Is  set  to  zero.  The  gun  tube  is  chosen  as  60  cm  long,  with  a 
projectile  offiiet  of  10  cm.  At  time  «  0,  the  gun  tube  Is  filled  with  hot  gas  at  a  pressure  of 
500  MPa.  The  resulting  charge  Is  762.8  g,  and  the  projectile  mass  was  set  at  5,000  grams 
(charge/mass  ratio  >  0.25). 

The  1-0  code  was  run  vfith  a  maximum  grid  spacing  of  5  cm,  2  cm,  and  1  cm.  Muzzle 
exit  takes  place  at  1.06  ms.  Figure  8  compares  the  velocity  profiles  some  time  after  muzzle 
exit  The  sharp  bend  In  velodty  near  the  tube  exit  Is  an  artifact  of  the  model.  In  the  model, 
steady-state  choked  flow  is  set  up  instantaneously.  This  leads  to  a  jump  in  velodty  from  the 
grid  point  just  inside  the  tube  to  the  grid  point  at  the  muzzle  exit  In  fact,  there  should  be 
some  delay  In  reaching  steady-state  conditions.  Otherwise,  agreement  Is  good,  aitiiough  not 
as  good  as  the  previous  case.  Rgure  9  compares  the  pressure  curves.  By  this  point  in  time, 
the  pressures  are  low.  Again,  agreement  Is  reasonable. 
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The  Comer  model  is  compared  to  the  1*0  simulation  with  a  grid  spacing  of  1  cm 
(Figures  10-15).  The  muzzle  velocity  and  average  gun  tube  temperature  at  muzzle  exit  are 
used  as  Input  into  the  Comer  model.  At  muzzle  exit,  the  1-0  code  does  show  a  linear  velocity 
profile.  The  pressure  profiles  are  ^jriy  flat  and  agreement  is  good.  Figure  12  compares  the 
velocity  profiles  some  time  after  muzzle  exit  The  comer  model  has  a  sharp  location  for  the 
raretection  wave,  while  the  1-0  code  shows  a  more  graduaJ  curvature.  The  agreement  near 
the  tube  entrance  Is  still  excellent,  and  the  location  of  the  rarefaction  wave  Is  about  correct 
The  pressure  curve  shows  similar  behavior.  At  a  later  time  (Rgures  14  and  15)  tire 
agreement  is  still  good.  So  the  blowdown  model  is  qualitatively  showing  the  proper  type  of 
behavior. 

5.  TEST  CASE 

The  thermal  management  code  was  first  applied  to  a  simulation  of  the  firing  of  Unicharge 
propellant  In  a  23-iitar,  52-call)er  cannon,  with  a  projectile  travel  of  6.9  m.  The  zone  six 
Unlchaige  produced  a  muzzle  velocity  of  about  945  m^  wHh  the  43.54  kg  (96  pound)  M549 
projectile.  To  make  a  fair  comparison,  a  liquid  propeHant  simulation  was  developed  that 
mimicked  as  doseiy  as  possible  the  Unicharge  simulation. 

As  a  starting  point  the  14>iiter  case  mentioned  in  Section  4.1  Is  used.  The  projectile  offset 
was  reduced  to  essentially  zero  (0.01  cm),  and  the  travel  extended  slightly  to  6.9  m.  The 
projectile  mass  was  decreased  siightiy  to  the  M549  projectile  mass.  The  projectile  resistance 
profiles  were  chartged  to  that  used  in  the  Urticharge  simulation.  Air  shock  in  front  of  the 
projectile  was  not  included.  Rnalty.  the  damper  was  actiusted  until  a  muzzle  velocity  near 
945  m/s  was  obtained. 

Again,  a  series  of  convergence  tests  were  made  (see  Table  2).  For  these  tests,  the 
integration  was  carried  out  to  50  ms.  Muzzle  exit  ocojrred  at  around  17.6  ms. 

Figures  16-19  compare  the  profiles  for  the  three  different  grid  sizes  (with  the  smaller  time 
step)  for  the  profiles  at  25  ms.  The  coarse  grid  Is  not  completely  accurate,  but  the  medium 
and  fine  grids  show  excellent  agreement 
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Rgure  12.  Simple  Blowdown.  Gun  Tube  Gas  Velodtv.  Tima  » 1.5  ms.  One-DImanslonal 
Model  (Une>.  Comer  Theory  (Pott. 


Figure  13.  Simple  Blowdown.  Gun  Tube  Gas  Pressure.  Time  =  1.5  ms.  One-DImenslonal 
Model  (Line).  Corner  Theory  (Dot). 
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Model  (Line^.  Comer  Theory  (Pott. 


Figure  15.  Simple  Blowdown.  Gun  Tube  Gas  Pressure.  Time  =  2.0  ms.  One-Dimensional 
Model  (Line).  Comer  Theory  (Dot). 
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Table  2.  Comparison  Between  Lumped  Parameter  and  1-D  Models  for  Thermal 
Management  Test  Case 


dxmax 

(cm) 

nmax 

Vm 

(m/s) 

Maximum 

Chamber 

Pressure 

(MPa) 

Run  time 
(hr/min/s) 

5 

50 

947 

273 

1:49 

5 

50 

945 

275 

3:42 

2 

200 

952 

274 

10:12 

2 

200 

947 

275 

27:53 

1 

500 

952 

275 

1K}9:09 

1 

500 

946 

276 

3:46:48 

.01 


.003 


.003 


.001 


.001 


03 


Grid  (Dot).  Fine  Grid  (Dash 
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Figure  19.  Test  Case.  Gun  Tube  Gas  Temperature.  Time  ■  25  ms.  Coarse  Grid  (Linel. 

Medium  Grid  rPott.  Rne  Grid  (Dashl. 

For  the  thermal  management  simulation  (Keller  et  al.  1991),  the  integration  was  carried  out 
to  500  ms.  This  was  to  make  certain  that  the  gas  in  the  gun  tube  cooled  to  near  ambient 
conditions.  Significant  heat  transfer  can  occur  long  after  muzzle  exit  as  long  as  the  gas  in  the 
tube  is  still  hot  For  this  simulation,  the  fine  grid  was  used  with  a  time  step  of  .0001  ms.  The 
grid  refinem^^nt  tests  reported  in  Table  2  indicate  that  this  is  actually  a  finer  grid  and  a  much 
smaller  time  step  than  required  for  numerical  convergence. 

Figure  20  shows  the  pressures  at  both  ends  of  the  gun  tube  as  a  function  of  time.  The 
pressure  essentially  reaches  atmospheric  pressure  shortly  after  100  ms.  Figure  21  shows  the 
corresponding  average  gas  temperature  in  the  gun  tube.  The  temperature  increases  as  the 
pressure  goes  up,  since  the  gas  is  being  compressed  and  hence  heated.  Once  the  pressure 
equilibrates  with  the  outside  air,  further  gas  temperature  decreases  are  very  gradual. 

Figure  22  shows  the  velocity  at  the  right-hand  end  of  the  tube  simulation.  This  is  the 
projectile  velocity  until  muzzle  exit  and  the  muzzle  gas  velocity  afterwards.  For  comparison, 
the  space  mean  sound  speed  in  the  tube  is  also  shown.  At  muzzle  exit,  by  coincidence,  the 
projectile  velocity  is  almost  the  same  as  the  local  speed  of  sound.  So  there  is  no  jump  in 
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300.0 


Rgure  20. 


Test  Case.  Pressura  at  Left  of  Tube  ^Line^. 
(Pott. 


Pressure  at  Right  of  Tu 


Figure  21.  Test  Case.  Average  Gun  Tube  Gas  Temperature. 
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Rgure22.  Teat  Case,  Gag  Velodtv  at  Right  of  the  Tube  fUne).  Averaao  Speed  of  Sound  ln 

velocity  as  stead^state  dwkecf  tk)w  is  establfaih^  The  flow  stays  ctiokeduntii  a  little  after 
100  ms.  When  the  flow  unchokes.  the  velocity  drops  more  rapidiy.  Because  of  inertia,  the 
pressure  in  the  tube  actuaUy  becomes  less  than  1  atm.  and  the  exit  velod^  becomes 
negative.  Gas  then  flows  in  and  out  of  the  tube  in  a  gradually  damped  fashion. 

Figure  23  shows  the  veiodty  profiles  in  the  tube  starting  shortly  after  muzzle  exit  The 
sharp  increase  in  veiodhr  just  at  the  odt  is  an  artifact  of  the  model,  since  the  exit  velodty 
instantaneously  reaches  ttw  steady-state  choked  flow  velodty.  The  rarefoction  wave 
propagates  very  siowiy  downtijbe.  since  the  local  gas  velocity  is  almost  at  the  speed  of  sound. 
Figure  24  shows  tire  velocity  profiles  at  more  widely  spaced  time  intervals.  The  rarefaction 
wave  is  damped  as  it  propagates  downstream.  This  is  tire  expected  behavior  (Comer  1950). 
Figures  25  and  26  show  the  corrr ''ponding  pressure  profiles. 
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6.  CONCLUSIONS 


A  regenerative  liquid  propellant  gun  code  which  is  1-D  in  the  barrel  and  lumped  parameter 
elsewhere  was  developed.  Before  muzzle  exit,  the  new  code  agrees  very  well  with  the 
entirely  lumped  parameter  model.  After  muzzle  exit,  the  results  are  qualitatively  as  expected. 
Combined  with  a  thermal  management  code  to  predict  wall  temperature,  the  model  allows  the 
detailed  simulation  of  barrel  heating  for  liquid  propellant  guns. 
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